Energy harvesting efficiency of piezoelectric flags in axial flows 

Sebastien Michelin^ 'B and Olivier Doare^'0 

^ LadHyX - Departement de Mecanique, Ecole polytechnique, 91128 Palaiseau Gedex, France. 
'^ENSTA Paristech, Unite de Mecanique, Chemin de la Huniere, 91761 Palaiseau, France 

(Dated: October 9, 2012) 

Self-sustained oscillations resulting from fluid-solid instabilities, such as the flutter of a flexible 
flag in axial flow, can be used to harvest energy if one is able to convert the solid energy into 
electricity. Here, this is achieved using piezoelectric patches attached to the surface of the flag 
that convert the solid deformation into an electric current powering purely resistive output circuits. 

■ Nonlinear numerical simulations in the slender-body limit, based on an explicit description of the 
coupling between the fluid-solid and electric systems, are used to determine the harvesting efficiency 

■ of the system, namely the fraction of the flow kinetic energy flux effectively used to power the 
^N) ■ output circuit, and its evolution with the system's parameters. The role of the tuning between 

the characteristic frequencies of the fluid-solid and electric systems is emphasized, as well as the 
O . critical impact of the piezoelectric coupling intensity. High fluid loading, classically associated 

■ with destabilization by damping, leads to greater energy harvesting, but with a weaker robustness 
to flow velocity fluctuations due to the sensitivity of the flapping mode selection. This suggests 

, that a control of this mode selection by a careful design of the output circuit could provide some 

opportunities of improvement for the efflciency and robustness of the energy harvesting process. 

V 

^ ■ I. INTRODUCTION 

C/3 1 The limited availability and environmental impact of fossile fuels motivate the development of renewable energy 
. ^ ■ sources. Significant research efforts are currently made to propose energy harvesting concepts and prototypes convert- 
c/2 ■ ing the kinetic energy of geophysical flows such as winds, rivers and oceanic or tidal currents into electricity [32j . In 
parallel, a particular attention is currently given to systems able to produce limited amount of energy from different 
vibration sources in order to power remote or isolated devices [29j . Classical fluid-solid couplings and instabilities 
such as vortex-induced vibrations, galloping and flutter in axial flows efltectively act as energy extraction mechanisms 
as they enable an energy transfer from the incoming flow to the solid body, and can therefore be used to produce 
_ electricity using displacement-based (e.g. electromagnetic converters) or deformation-based (e.g. piezoelectric mate- 
^ . rials) conversion mechanisms [SS US, 113] ■ Because they are based on fundamentally-different mechanisms, such flow 
energy harvesters may be attractive complements to the existing wind- and water-turbines technologies, and properly 
assessing fundamental upper bounds on their respective efficiency is therefore of critical importance. 



^vq ■ A flexible plate placed in an axial flow becomes unstable to flutter above a critical flow velocity when the destabilizing 
—-^ i pressure forces dominate the stabilizing effect of the structure's rigidity [l^llSH^- This critical velocity depends on 

■ the plate's properties (e.g. density, size and rigidity) and can therefore be adjusted in the system's design to be lower 
] than the typical flow velocity. This so-called flapping flag instability leads to self-sustained large-amplitude flapping 

of the plate in the form of traveling bending waves J22[|. that ca n be used to produce electricity using, for 

example, piezoelectric patches attached to the plate's surface [ul^. [isl [T7| . 

An important research effort is required in order to assess the amount of energy that can be harvested using such 

■ devices and investigate possible intrinsic limits or potential optimization strategies of their efficiency. In a theoretical 
H . or numerical framework, the conversion mechanism and output circuit must be described, to properly include the 

■ coupling of the fluid-solid and electric systems. Energy harvesting eventually amounts to an extraction of energy from 
the solid dynamics. Hence, a first and simpler model for the harvesting mechanism is an additional structural damping 
(e.g. Kelvin- Voigt), and assessing the s yste m's efficiency is then equivalent to determining how much energy can be 
dissipated by the flapping structure l27l. [30| . Indeed, increasing damping would lead to a larger energy dissipation 
but eventually will re-stabilize the system and reduce its harvesting efficiency. Although simple to implement, this 
representation is not complete as it assumes that energy is instantaneously and immediately dissipated and can not 
represent the dynamics of the electrical circuit or of the coupling mechanism. 
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The originality of the present work is to propose instead a fully-coupled description of a fluid-solid-electric system, 
namely a flexible plate in axial flow covered with piezoelectric patch pairs powering simple resistive elements. Recently, 
Doare & Michelin [l^l followed this approach to study the impact of the piezoelectric coupling on the linear stability 
of a two-dimensional plate and on the solid-electric energy transfers. In particular, the role of the tuning of the 
fluid-solid and electric characteristic time-scales was emphasized, and a destabilization by the piezoelectric coupling 
was identified in the case of large fluid loading, associated with the destabilization by damping of negative energy 
waves 0, [ll[ . The present study extends this approach to study numerically the nonlinear dynamics of this fluid-solid- 
electric system in the case of a slender flexible plate, and to determine its harvesting efficiency. Here, the system's 
efficiency is defined following the classical definition used for wind-turbines, as the ratio of the mean power output 
and of the mean kinetic energy flux through the section occupie d by the device in the flow. In that sense, it differs 
from the measures of efficiency used in other existing studies [l^ [301 ■ 

In f|ni the model used to describe the dynamics of the flapping piezoelectric flag is presented. Section Hill presents 
a short summary of the linear stability results in the case of a slender plate. In mVl the numerical solution of the 
coupled dynamics is addressed, the system's efficiency is defined and the impact of the different system parameters 
on this efficiency is discussed. Finally, conclusions and perspectives are presented in jVl 



II. PRESENTATION OF THE FULLY-COUPLED MODEL 



A. Piezoelectric flag dynamics 



The system considered here consists of a rectangular flexible plate of length L, width H and thickness h {h <C H, L) 
placed in a steady flow of density p and velocity Uoo- The plate is inextcnsiblc and clamped at its leading edge; for 
simplicity, only purely planar motions of the plate are considered, so that the plate's position X is only a function 
of the streamwise curvilinear coordinate S and time T, and the solid does not experience any spanwise displacement 
nor twist. The local orientation of the flag with respect to the horizontal axis is noted 9{S,T) (Figure [1]). In the 
following, lineic quantities will be defined per unit length in S. 

The surface of the plate is covered by pairs of piezoelectric patches (Figure [Ud) with streamwise length I <C L 
and width H. The negative electrodes of each patch are shunted through the plate and the positive electrodes are 
connected to the output circuit. The deformation of the flag is coupled to the output circuit through the piezoelectric 
coupling: (i) stretching and compression of the patches due to the local curvature induces charge transfers between 
each patch's electrodes and (ii) an electric voltage applied to its electrodes results in an additional internal torque 
on the piezoelectric patch and on the flag. Considering the limit of a continuous coverage by patches of infinitesimal 
length [1, , the local electric state can be described in terms of the electric voltage between the positive electrodes 
of each patch, V{S,T), and the charge transfer Q{S,T) per unit length in the streamwise direction. In this limit, 
which differs from the single-patch approach of Akcabay & Young [l[ , both quantities are continuous functions of S 
and T, and the piezoelectric coupling imposes that: 

an 

Q = CV + X*gg, (1) 

M = B—-x*V, (2) 

where A4 is the total internal torque in the piezoelectric flag, and c and x* ^-rc the lineic capacitance and piezoelectric 
coupling coefficient, directly related to the material and geometric properties of the patch pair [l^ . An Euler-BernouUi 
model is assumed for the dynamics of the piezoelectric flag with B the effective flexural rigidity of the three-layer 
piezoelectric plate [for more details, see [l^ Il9| . 

The positive electrodes are connected to a purely resistive circuit of lineic conductivity g (Figure [2^), such that 

11+5^^ = 0. (3) 



The conservation of momentum and inextensibility condition for the flag lead to: 

Ffluid, (4) 
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FIG. 1: (a) Slender flexible plate flapping in a uniform axial flow, (b) Two-dimensional flapping of a slender flexible plate 
covered with pairs of piezoelectric patches. 




FIG. 2: (a) Piezoelectric patch pair powering a purely resistive circuit, (b) The piezoelectric patch pair is characterized by the 
current generated by its deformation and by its capacitance c. 



with /i the lineic mass of the piezoelectric flag, Ft{S, T) the local tension, acting as a Lagrangian multiplier to enforce 
the plate's inextensibility (O, and the internal piezo-elastic torque in The following clamped- free boundary 
conditions must also be satisfied; 



M 



X = 0, 
dM 



as 



61 = 
Ft = 



at 5 = 0, 
at S* = L. 



The conservation of mechanical and electrical energy takes the following form: 



where 



5fc = 



r 



dT 



9X 



dT 



dT 



dS, 



S'e.l = 



-cV'^dS 
2 



(6) 
(7) 

(8) 

(9) 



are respectively the kinetic and potential elastic energy of the flag, and the energy stored in the capacitance of the 
piezoelectric elements, and 



fluid • 



-X 



(10) 



are the rate of work of the fluid forces, the rate of energy transfer from the solid to the electric circuit and the power 
used in the output circuit, respectively. For a purely resistive circuit ^ is always strictly positive, and in permanent 
periodic regime, (^) = (^) = (/^), with (•) the time-averaging operator. 
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B. Fluid modeling - Lighthill's theory 

The relative motion of the solid body with respect to the incoming flow results in fluid forces Ffiuid applied on 
its surface. In the particular limit of a slender body (H <C L) and for a purely potential flow, the extension of 
Lighthill's Elongated Body Theory to large amplitude displacements leads to the following leading order expression 
for the reactive fluid forces Freac associated with the local transverse motion of each cross section along the plate: 

with Tfia the non-dimensional added mass coefficient of the local cross-section, namely th^ — tt/A for a flat plate. In 
((TT|) , Ut and Un are respectively the tangential and normal components of the local relative velocity of the solid with 
respect to the incoming flow: 

U = — - C/ooe:r = C/re^ + C/„e„. (12) 
ol 

Initially proposed by Lighthill [2l|, this so-called Large Amplitude Elongated Body Theory (LAEBT) was recently 
shown to provide a good estimate of the transverse fluid forces, in comparison with RANS simulations on a towed 
and deforming fish body 0]. However, Candelier et al. Q emphasized that this purely reactive formulation can not 
by itself represent properly the deformation amplitude of freely-moving bodies, as such effects as drag and separation 
will be significant and must be accounted for by an additional resistive component Fj-osist [see for example |3l| . In 
the case of a freely-flapping slender body, Singh, Michelin & de Langre [1^ indeed observed that the purely reactive 
model would lead to non-physical overestimates of the flapping amplitude. Following Eloy et al. [l^ and Singh et al. 
[13], the present model only retains the resistive drag associated with the plate's normal displacement 

Frcsist = -^pHCoUn |C/„| e„, (13) 

with Cd = 1.8 for a flat plate in transverse flows. 

The reactive part of the LAEBT corresponds to the asymptotic limit of the potential flow equations when H/L <C 1 
0, but the recent work of Eloy et al. ^ showed nonetheless, using comparisons with wind-tunnel experiments, that 
the combination of the reactive and resistive components ()lip and (|13|) can provide a good prediction of the flapping 
properties of the plate even when H/L = 0{1). In the following, an aspect ratio H* = H/L = 0.5 will therefore be 
considered. 



C. Non-dimensional equations 



Equations p])- P^ are non-dimensionalized using L, L/Uoc, pHL^, f/ooi/M/c and Uod^/W^ as characteristic length, 
time, mass, voltage and charge density, respectively: 





d 


dt^ ~ 


ds 


/fluid = 
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Jrer - TT 



1 de 
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ds \ U*^ ds U 
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and the tension /^ is obtained using the inextensibility condition [see for example [H, [l^l 

9x 

ds 

The clamped-free boundary conditions become 

at s = 0, X = 0, 6* = 

at s = 1, fr = TT - ctU V ^ — ~ aU — = 0. 

ds ds'' ds 



(14) 

(15) 
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Five non-dimensional parameters characterize the system, namely the fluid-solid inertia ratio, the non-dimensional 
velocity U* , the coupling coefficient a, the tuning coefficient of the fluid-solid and electric system /3 and the aspect 
ratio of the plate H* : 

The originality of the present work is to offer a full description of the fluid-solid-electric system. Equations (|14p , 
and (|17p show that the effect of the piezoelectric coupling is more complex than the simple Kelvin- Voigt damping 
model generally assumed for simplicity in most studies on energy harvesting flags [2^, |30|. Indeed, combining 
and (|17p. one obtains 

Bv + v^ ^TT^. 22) 

Equation (22) shows that the effective damping introduced by the piezoelectric is frequency-dependent. In fact, a 
Kelvin- Voigt damping model could only be recovered in the particular limit of /3 <C 1 and finite aP/U* . However, this 
asymptotic limit is unlikely to be achieved in practice because of the material restrictions on the coupling coefficient 
a for currently-available piezoelectric materials [T2j . 

III. LINEAR STABILITY ANALYSIS 

The linear stability of the piezoelectric flag is first analyzed to identify the impact of the piezoelectric coupling 
and output circuit on the stability properties of the system, and also identify the operating regime of the harvesting 
devices, namely the parameter values for which self-sustained oscillations can develop. The present linear study only 
differs from that in Doare & Michelin [l^l by the fluid model considered, that corresponds to a different range for 
the plate's aspect ratio, therefore only the main results will be reminded and the reader is referred to this previous 
contribution for more in-depth analysis of the linear stability. 

The displacement of the flag is purely vertical and noted y{s,t) ^ 1. At leading order, P^ - ipUl) simplify into the 
following linear systems for (y,w): 

(1 + ™.M-H-)|i + 2™«M-ff-^ + ™.M'//-g| + J^g - = 0, (23) 



with boundary conditions: 



dt U* ds^dt 



ats = 0, y=^ = (25) 



s 



ats = l, -Ji-aU*v^-^-aU*—=0. 26 
os^ os'^ OS 

Searching for solutions of the form [y,v] ~ Re F]e~"^*^ , become an eigenvalue problem for lu and 

[Y{s),V{s)], that is solved numerically using a Chebyshev collocation method to determine the stability of the 
piezoelectric flag, and in particular the critical velocity above which the flag becomes unstable (Figure [3]). 

The piezoelectric coupling a enables the transfer of energy from the fluid-solid system to the electrical circuit 
where part of it is dissipated, resulting in a net damping on the solid motion. This additional damping is therefore 
expected to increase the critical velocity in comparison with the uncoupled flag (a = 0), an effect indeed observed for 
^* ^ 1 (Figure 13^). At larger M*, the piezoelectric coupling instead destabilizes the system, at least initially. This 
destabilization by damping was previously reported in the case of a two-dimensional flag by Doare & Michelin p^ . 
and is associated with the existence of negative energy waves in the local stability analysis of the non-dissipative flag 
. From an energy harvesting point of view, it increases the operating range of the piezoelectric flag as self-sustained 
oscillations develop for lower velocities. 

For a flxed piezoelectric coupling, f3 measures the tuning of the fluid-solid and electric time-scales of the system. 
When forced by the flag at a frequency much lower than 1/trc — g/c (/? ^ 1), the output resistive elements are 
seen by the piezoelectric as short circuits, and the voltage at the electrodes remains negligible. The critical velocity 
is therefore equal to that of the uncoupled piezoelectric flag (a = 0) as no piezoelectric feedback is applied on the 
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FIG. 3: (Left) Critical velocity threshold U* as a function of the mass ratio M* for (3 = 1 and varying a. (Right) Critical 
velocity threshold as a function of the mass ratio M* for a = 0.5 and varying j3. 



structure. For a large forcing frequency (/3 3> 1) however, the resistive elements are seen as open circuits, and from 
Eq. ([1]), the voltage at the piezoelectric's electrodes is proportional and opposite to the curvature: the piezoelectric 
coverage then acts as an additional rigidity on the system. Between these two limit regimes, a destabilization is 
observed for large M* which corresponds to the destabilization mechanism mentioned above (Figure [3]3) . 

These results confirm and extend to the slender-body limit the conclusions of the infinite-span flag analysis of Doare 
& Michelin [l^ . It is worth noting that the results obtained with both models differ mostly at low M* , consistently 
with the results of Eloy et al. [l^ on the impact of aspect ratio on the flag stability. 



IV. NON-LINEAR DYNAMICS OF A PIEZOELECTRIC FLAG 



To determine the amount of energy that can be produced using such a system, the nonlinear dynamics of the 
piezoelectric flag must be studied, in particular to determine its flapping amplitude and frequency. 



A. Non-linear simulations and energy harvesting efficiency 

Following Alben the nonlinear system (jl4p - (|20p is integrated numerically in time using a second-order accurate 
implicit method, and spatial derivatives are computed using Chebyshev collocation. Starting from rest {9{s,t < 0) = 
0), the flag is excited by a small perturbation in the vertical component of the upstream flow. The harvested energy 
is computed as the temporal average of the non-dimensional power V — ^/{pU^HL) dissipated in the resistive 
elements in permanent regime: 

fn the previous equation, the temporal average is understood and computed as follows: when the system converges to 
limit-cycle oscillations, it is defined as the mean value over a period of oscillation, but when no limit-cycle oscillation 
can be identified, it is computed as the statistical average of ^ over a long enough time frame. Similarly, the non- 
dimensional flapping amplitude A = / L is defined from the trailing edge displacement ye(t) as a measure of the 
peak flapping amplitude: 



A=Vm)- (28) 

The harvesting efficiency of the system, rj, is defined as the fraction of the fluid kinetic energy flux through the 
cross-section 2s/ H occupied by the flag (Figure [T]) actually transferred to the output circuit, namely 



yUl, X 2si^H ~ ~A' ^^^^ 
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B. Non-linear flapping dynamics 

Above the critical velocity U*, defined using the linear stability analysis of Section IlIIl an initial perturbation 
of the flag's state of rest leads to an exponential growth of the flapping amplitude until saturation is reached, and 
the permanent regime takes one of the two following forms: (i) a strongly periodic regime characterized by the 
identification of a limit-cycle in phase-space or (ii) a more complex non-linear regime where no clear limit-cycle can 
be identified. This transition from periodic to non-periodic regime has been observed in numerous experimental [T^ 
and numerical studies [3, liol . [22j . and has been conjectured to result from the non- linear interactions of different 
fundamental modes. Limit-cycle oscillations are particularly interesting from an energy harvesting point of view as it 
provides steady output current amplitude and frequency. 

Even below the transition to chaotic flapping, non-periodic flapping regimes can be observed as the systems switches 
from one flapping mode to another when one of the parameters (e.g. U*) is modified. This mode switch results in a 
change of flapping amplitude and frequency, but also of the flag kinematics resulting in a modification of the forcing 
distribution on the piezoelectric elements (Figure |4]) . 

C. Variations of the energy harvesting efficiency 

In the following, the impact of the different parameters on the harvesting efficiency is presented. 

1. Effect of the tuning ratio 

The ratio /3 = cUoo/{gL) measures the relative tuning of the ffuid-solid and electric time scales, Tadv = L/Uoo and 
Tfic = c/g, respectively. All other parameters being fixed, it is observed that the harvested energy efficiency reaches 
a maximum when puj = 0(1) where lo is the non-dimensional flapping frequency of the flag (Figure [S]). 

The existence of this maximum comes as no surprise: when /3 <C 1 and /? ^ 1, the resistive element acts as a 
short-circuit or open-circuit, respectively. In both cases, no energy is dissipated and 77 = 0. /3a; = 0(1) corresponds 
to a forcing of the RC output circuit at its characteristic time-scale, which is expected to result in maximum energy 
dissipation in the resistive element. The forcing frequency is however not a property of the ffuid-solid system only, 
but is instead the result of the nonlinear coupling between the ffuid-solid system and the electric output through 
the piezoelectric material (Figure [5]). Similarly, the flapping amplitude is signiflcantly modified when (3 is varied: in 
particular, for M* = 1 (Figure [5^), a sharp drop in the flapping amplitude is observed as /3 is increased. 

This result is also confirmed on Figures IH] and [T] For each value of U* , an optimal tuning ratio can be determined 
and the optimal /? is a decreasing function of U* (Figure [6)3). This is consistent with the observed increase in flapping 
frequency w with U* (Figure [7]) and the criterion fjoj = 0(1) for optimal energy harvesting. 

2. Effect of the flow velocity 

Previous experimental results on the dynamics of ffexiblc flags have established that the flapping amplitude is in 
general an increasing function of the non-dimensional velocity U* above the instability threshold [l4l . [iH l25l | , before 
saturation of this flapping amplitude is reached. For a given flapping mode shape and frequency, the harvested power 
V varies quadratically with the amplitude A, therefore it is expected that raising U* will lead to an increase in the 
system's efficiency. This is confirmed partially in Figure [SI when a given flapping mode remains dominant, 77 is indeed 
an increasing function of U* , mainly due to the associated increase in flapping amplitude. However, when a mode 
switching event occurs as described in Section llVBl a sudden decrease of the efficiency is observed, mainly associated 
with a reduction in the flapping frequency (Figure [7]). 

Figure [7] shows that the nonlinear flapping frequency is very close to the frequency of one of the unstable linear 
modes of the piezoelectric flag. A mode switching event, as U* is increased, consists of a transition from one linearly 
unstable mode to another with lower frequency. A study of the associated linear growth rate however docs not show 
any coincidence of such event with a change in the most unstable linear mode, and this mode switching event is 
therefore the result of a purely nonlinear mechanism. Figure [5] shows that such mode switching events take place at 
lower values of U* for lighter flags (large M*) while for M* < 1, no such even is detected below U* = 20. 

This mode selection mechanism is also observed for a flapping flag without any piezoelectric (a = 0). Regardless 
of its origin, its importance is however essential for the performance of the energy harvester: as long as the same 
nonlinear flapping mode can be maintained, the efficiency of the system increases with U* and the occurrence of 
a mode switching event results in an important performance loss for the device. A better understanding of this 
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(a) [/* = 10.5, 77 =3.8% (h) U* = U, rj =1.7% (c) t/* = 14.5, t? =3.3% 




FIG. 4: Limit-cycle oscillations: (Top) Flapping mode shape, (Center) Time-series of the non-dimensional harvested power 
'P(f) and (Bottom) Phase-space trajectory for the trailing edge orientation 6e{t) for M* = 10, a — 0.5, /3 = 0.158 and (a) 
U* = 10.5, (b) 11 and (c) 14.5 from left to right. 



phenomenon and, in particular, of the impact of the piezoelectric coupling on the transitions, is therefore required 
and could lead to significant improvements of the harvesting efficiency by constraining the system to a more efficient 
flapping. 

3. Effect of the mass ratio 

The linear analysis of Doare & Michelin [l^ identified significant differences in the performance of lighter (large 
M*) or heavier fiags (small M*), as measured by the energy transfer from the structure to the output circuit. Higher 
performance at large M* was associated with the destabilization by damping of negative energy waves. 

A similar result is observed here in nonlinear simulations for the harvesting efficiency rj (Figures |8] and [9]): harvesting 
efficiencies up to 10-12% can be achieved for M* ~ 20 and U* < 20, while the optimal value of rj is less than 1% 
below M* = 1. 

Comparing Figures |Sl[a) and (b), the optimal M* appears as closely related to the tuning parameter /3, emphasizing 
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(a) U* = 14.5, M* = 1 



(b) U* = 14.5, M* = 10 



FIG. 5: Evolution with the tuning ratio /? of (Top) the non-dimensional flapping frequency uj and (Bottom) the non-dimensional 
flapping amplitude A (dashed), harvested energy Q (dash-dotted) and harvesting efficiency rj (solid) for a = 0.5, U* = 14.5 
and (a) A'P = 1 or (b) Af* = 10. On the top, the dotted line correspond to the variations of 



again the importance of the synchronization of the fluid-solid and electric systems: for small /S, regions of greater M* 
will be optimal as they correspond to larger flapping frequencies while heavier flags (smaller M*) will be optimal for 
larger values of /3. 

Finally, Figure [5] shows the optimal-tuning efficiency as a function of (Af *,[/*). Up to 12% of the kinetic energy 
flux can be harvested for the largest value of M* and U* considered. However, it is also important to emphasize that 
this parameter region corresponds to closely-spaced mode switching events, making the efficiency of the system quite 
sensitive to fluctuations in the flow velocity. 



4- Effect of the piezoelectric coupling 



The coupling coefficient a is a measure of the intensity of the fluid-solid and electric systems' forcing on each other, 
and as such is clearly expected to impact the amount of energy transferred to the output load. Figure [TU] shows 
the evolution of A and rj when a is increased. For small coupling a <g; 1, the flapping dynamics is only marginally 
modified and the amplitude of the charge transfer q and electric potential v increase linearly with a as seen in Eq. (|16p . 
As a result, Q and rj initially increase quadratically with a (see inset on Figure [TU)) . However, when a is increased 
further the feedback piezoelectric coupling modifies the flapping dynamics resulting in a linear decrease of the flapping 
amplitude and harvesting efficiency and, eventually, the restabilization of the system. One can therefore identify an 
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10"^ 10"^ 10° 10^ 10"^ 10"^ 10° 10^ 

/3 /? 

(a) M* = 1 (b) M* = 10 

FIG. 6: (Top) Harvesting efRciency rj and (Bottom) flapping amplitude ^ as a function of the tuning ratio /3 and the non- 
dimensional flow velocity U* for (a) M* = 1 and (b) M* = 10. For both cases, a = 0.5. The black dashed line corresponds to 
the instability threshold [/* below which rj = A = Q. 



optimal value of the coupling coefHcient, in the same way that an optimal damping was determined for maximum 
energy dissipation in Singh et al. [l^, HI]. The value of the optimal coefficient Uc clearly depends on the other 
system parameters, and will be greater when the flag is far from its stability threshold before piezoelectric coupling 
is introduced, or when destabilization by damping occurs as for larger M* . 

Achieving the optimal a is however not necessarily possible practically: a is a characteristic of the material's 
electric and mechanical properties and is of the order a ~ 0.3 for typical piezoelectric materials such as PZT and even 
lower for PVDF [l2|. Except in the vicinity of the instability threshold, the optimal a leading to maximum energy 
efRciency is however expected to be greater than this value, suggesting that an optimization of the piezoelectric flag 
design or future technical improvements in the properties of available piezoelectric materials can potentially increase 
the achievable values of a and lead to significant efficiency gains. 

V. DISCUSSION AND PERSPECTIVES 

The present study focused on the fully-coupled dynamics of a classical fluid-solid system, a flexible plate in axial flow, 
and a simple resistive circuit coupled through piezoelectric patches attached to the surface of the plate and converting 
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5 10 15 20 5 10 15 20 

u* u* 

(a)/3 = 0.03 (b)/3 = 0.31 



FIG. 7: Evolution with U* of the limit-cycle dominant frequency ui (black star) for M* = 10, a = 0.5 and (a) /3 = 0.03 and (b) 
(3 — 0.31. In each case, the frequencies of the different linear modes are shown: light dotted lines correspond to stable modes 
and thick grey lines to unstable modes. On each figure, from left to right, unstable frequencies correspond to flapping modes 
of increasing order and decreasing characteristic wavelength. 




FIG. 8: Harvesting efficiency 7] as a. function of the mass ratio M* and normalized velocity U* for (3 — 0.1 (left) and /3 = 1 
(right). In both cases, a = 0.5. The black dashed line corresponds to the instability threshold U* below which r; = 0. 



the plate's bending deformation into an electric current. In the limit of continuous coverage by infinitesimal patches, 
the energy harvesting efficiency was determined as a function of the different system parameters, namely the inertia 
ratio, the non-dimensional flow velocity, the coupling coefficient and the tuning ratio. For realistic coupling coefficients, 
as much as 10% of the kinetic energy flux can be transmitted to the output circuit, but this efficiency was found to 
be highly sensitive to several important parameters, in particular the coupling coefficient and the flow velocity. 

This study confirms the results by Doare & Michelin [l2| on the impact of destabilization by damping and on the 
solid-electric energy transfers: in nonlinear saturated regimes, those parameter regions correspond indeed to maximum 
energy harvesting efficiency. The critical role played by the tuning ratio is also confirmed: maximum energy transfers 
are obtained when the output circuit characteristic timescale is tuned to the flapping frequency. This frequency is 
itself determined through the nonlinear coupling of the fluid, solid and electric systems, and modifications in the 
fiapping frequency associated with a switch in the flapping mode directly impact the efficiency of the system and 
its robustness to fluctuations in flow velocity for example. Controlling the flapping mode selection is therefore an 
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FIG. 9: Harvesting efficiency obtained for tfie optimal tuning /3 as a function of the mass ratio M* and non-dimensional velocity 
U* for a — 0.5. For each value of M* and (7*, /3 is chosen so as to maximize the harvesting efficiency. The black dashed line 
shows the minimum of the instability threshold over all possible values of P, below which 77 = for all /3. 




FIG. 10: Evolution of the flapping amplitude A (dashed) and harvesting efficiency rj (solid) with the coupling coefficient a for 
P — 1, M* — 0.5 and U* = 15. The inset show the efficiency's scaling with a in the limit of a <C 1 and Qc — Q <?; 1 with 
Qfc ~ 0.82 the critical value of a leading to restabilization of the piezoelectric flag for those parameter values. 



important challenge for the improvement of the efficiency of this model energy harvester, and should be considered 
in future work. 

By coupling the fluid, solid and electric systems in a nonlinear model and by using an explicit description of the 
coupling mechanism and output circuit, the present approach provides some important insight on the nature and 
importance of the feedback of energy harvesting on the solid dynamics, as illustrated by the modification of the 
flapping amplitude and frequency, for example. Even though the simplest possible circuit (a purely resistive element) 
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was used here, the impact of the tuning ratio /3 on the efficiency suggests that significant efficiency gains should be 
expected through a careful design of the output circuit, using more complex and possibly active circuits, as well as 
state-of-the-art power electronics techniques such as synchronized switching techniques |20l |. 
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